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ABSTRACT 

Aims. We present colour-magnitude diagrams (CMDs) of young massive star clusters in several galaxies located well beyond the 
Local Group. The richness of these clusters allows us to obtain large samples of post-main sequence stars and test how well the 
observed CMDs are reproduced by canonical stellar isochrones. 

Methods. We use imaging of seven clusters in the galaxies NGC 1313, NGC 1569, NGC 1705, NGC 5236 and NGC 7793 obtained 
with the Advanced Camera for Surveys on board the Hubble Space Telescope and carry out PSF-fitting photometry of individual stars 
in the clusters. The clusters have ages in the range ~ (5 - 50) x 10^ years and masses of ~ 10'Mq-10*Mq. Although crowding prevents 
us from obtaining photometry in the inner regions of the clusters, we are still able to measure up to 30-100 supergiant stars in each 
of the richest clusters. The resulting CMDs and luminosity functions are compared with photometry of artificially generated clusters, 
designed to reproduce the photometric errors and completeness as realistically as possible. 

Results. In agreement with previous studies, our CMDs show no clear gap between the H-burning main sequence and the He- 
burning supergiant stars, contrary to predictions by common stellar isochrones. In general, the isochrones also fail to match the 
observed number ratios of red-to-blue supergiant stars, although the difficulty of separating blue supergiants from the main sequence 
complicates this comparison. In several cases we observe a large spread (1-2 mag) in the luminosities of the supergiant stars that 
cannot be accounted for by observational errors. We find that this spread can be reproduced by including an age spread of ~ (10 - 
30) X 10* years in the models. However, age spreads cannot fully account for the observed morphology of the CMDs and other 
processes, such as the evolution of interacting binary stars, may also play a role. 

Conclusions. Colour-magnitude diagrams can be successfully obtained for massive star clusters out to distances of at least 4-5 Mpc. 
Comparing such CMDs with models based on canonical isochrones we find several areas of disagreement. One interesting possibility 
is that an age spread of up to ~ 30 Myr may be present in some clusters. The data presented here may provide useful constraints on 
models for single and/or binary stellar evolution. 

Key words. Open clusters and associations, galaxies: star clusters, galaxies: spiral, stars: Hertzsprung-Russell and C-M diagrams 



. !^ ' 1 . Introduction leal evolution) it is therefore essential to strive towards the best 

possible understanding of the properties of individual stars. This 
^ Stars profoundly affect the properties of the galaxies they m- ^^^^ ^^jy integrated light measurements, but also for 
d habit. They are the primary sites of nucleosynthesis and drivers i^odelling of colo ur-magnitude diagrams (CMDs) of resolved 
of chemical evolution and they conti'ibute to the energy budget ^^^u^ populations dGallart et aljr2 005). 
of the interstellar medium through feedback processes such as , , ^ , . . . , 

supernovae and stellar winds. In most galaxies they also deUver , Although stars spend most of their lifetimes in the core- 
most of the ultraviolet, optical and near-infrared light. In order ^ydrogen burning phase on the main sequence (MS), post- 
to correctly interpret observations of galaxies and constrain their ^S stai-s provide a significant fraction of the bolometnc lumi- 
past evolution (e.g. in terms of star formation histories and chem- ""^i^y °f f ^^^"^ population at all except the very youngest 

(few Myr) ages (e.g. Renzini & Buzzoni 1986). Intermediate- 
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* Based on observations made with the NASA/ESA Hubble Space ^uch stars is particularly uncertain. One example is the relative 

Telescope, obtained from the data archive at the Space Telescope numbers of blue and red core-He burning stars, which is very 

Institute. STScI is operated by the association of Universities for sensitive to assum ptions about input pai-ameters such as mass 
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the red/b lue ratio e.g. as a functio n of metallicity seen in star 
clusters ("Eggenbe rger et all l2002h . although somewhat better 
agreement is fo und when comparing trends in this rat io ver- 
sus a ge (fPohm-Palmer et aP Il997t Dohm-Palmer & SkillmanI 
Another issue concerns the extent of the "blue loops" in 
the core-He burning phase, which depends on the detailed struc- 
ture and chemical profile of a star These, in turn, are affected 
e.g. by the treatment of convective overshooting (Matraka et al. 
[1982 ; Ritossa|11996|), with the net effect that more overshooting 
tends to lead to less extended (redder) blue loops. For standard 
assumptions, models tend to underpredict the extent of the blue 
loops. While most stellar evolution models predict a gap in the 
colour-magnitude diagram (CMD) between the H and He burn- 
ing stars (the "Blue Hertzsprung Gap" or BHG; IChiosil [19981 
ISalasnich et al.|[T999l) . the post-MS radius evolution is very un- 
certain for massive stars and star cluster C MDs often show sig - 
nificant numbers of stars in the BHG (e.g. [Fischer etan[T99l . 
The nature of these stars is poorly understood an d it is un- 
clear whether they are even core H or He-burning (IVink et al.l 
although in some cases their CNO abundances suggest 
that they have underg one a first dredge-up as red supergiants 
(ISalasnich et alJ[T999l) . Some models do predict more extended 
blue loops, but this requires that a very large fract ion of the H- 
rich envelope is lost in the red supergiant stage (iHirschi et al.l 
120041: ISalasnich et al.l[T999l) . This may not be compatible with 
the fact that SN 1987 A still had a massive H-rich envelope when 
it exploded as a blue supergiant ( lWooslevlll988l) . Complicating 
matters further a large fraction of massive stars are m embers of 
binary systems jSana et al. [ [2008L[\te3?7^td][2009h . many of 
which may interact or even merge during a phase of common en- 
velope evolution (T aam & Sandquist 2000; Eldridge et aj. 2008; 
ISansom et al ] [2009h . 

Historically, star clusters have long served as important 
testbeds for stellar models. A prime reason for this has been 
the assumption that all stars in a cluster are very nearly co- 
eval, so that a cluster CMD can be directly compared against a 
model isochrone. Furthermore, under this assumption all post- 
MS stars in a given cluster will have very similar (initial) 
masses, so that the step of interpolation in the tracks to pro- 
duce isochrones can often be eliminated for compariso n with ob- 
serve d CMDs or calculation of integrated properties dMarastonl 
[T998h . Models of low-mass stars in particular have been ex- 
tensively tested by comparison with globular cluster CMDs 
( Ren zini & Fusi Pecci 1988; Palmieri et al. 2002; Salaris et al. 
r2007h . in some cases providing evidence of multiple stellar pop- 
ulations (|Piottoll2009l) . Similarly, high precision photometry has 
shown that the CMDs of several intermediate-age (~ 1 - 2 Gyr) 
clusters in the Magellanic Clouds are not well fitted by standard 
isochrones, but instead show a significant spread in the main se- 
quence turn-off magnitudes and luminosities of red clump stars. 
One possible explanation f or these fe atures is a significant age 
spread in these clusters (Mac kev et al. l 2008; Milone et al. 2009), 
although other effects (such as stellar rotation) may also play a 
role (Bastian & de Mink 2009; but see also Gh'ardi et al.ll201 ll) . 
However, interpretation of cluster CMDs remains easier than 
the complex mix of stellar populations that make up the field, 
and clusters remain attractive targets for comparison with stellar 
models. In fact, such a comparison may yield important infor- 
mation on possible age spreads within clusters, in addition to 
providing constraints on the stellar models. 

One difficulty is that a typical low-mass young (open) clus- 
ter contains no or only a few post-MS stars. As an illustration. 
Fig. [T]shows the number of stars with \ogg < 3.0 vs. age, pre- 
dicted for a single-age stellar population with an initial mass of 



10^ Mq. We have assumed a Salpeter IMF over th e mass range 
0.15 < M/Mq < 100 and used Padua isochrones (Bertelli et alj 
2009) to obtain \ogg values versus age and (initial) stellar mass. 
The log g limit a dopted in this figure co mfortably includes all lo- 
cal supergiants dLvubimkov et al.|[2010i) . The figure shows that, 
for ages less than about 30 Myr, only 20-30 such stars are ex- 
pected per 10^ Mq, with an increase at older ages. The figure 
is consistent with observations of post-MS stars in massive star 
clusters, e.g. of the rich LMC cluster NGC 1866 which has 100 
post-MS stars, an age of about 160 Myr and a mass o f about 
1.3 X lO^Mo ([Fischer et alj [19921; iBarmina et al.ll2002l) . or the 
"Red Supergiant Cluster 2" (age ~ 10 Myr, initi al mass ~ 400 00 
Mq) which hosts 26 red supergiant stars (Dav ies et al.[[2007[) . It 
is thus clear that cluster masses of at least several times W^Mq 
are required in order to get significant samples of massive post- 
MS stars. Such clusters are rare in the Milky Way and even in 
the Local Group, making it difficult to obtain samples that cover 
a range in age (or, equivalently, main sequence turn-off mass), 
metallicity, etc. 

If the search volume is extended beyond the boundaries of 
the Local Group several galaxies hosting young clusters with 
M <: 10^ Mq become available. The trade-off is that photometry 
of individual stars becomes extremely challenging due to severe 
crowding, even if imaging from the Hubble Space Telescope 
(HST) is used. One of the closest examples of a galaxy host- 
ing s uch massive, young star clusters is NGC 1569 at ~ 3.4 
Mpc ([Grocholski e t al. 2008). Using imaging obtained with the 
High Resolution Channel (HR C) of the Advanced Camera for 
Surveys (ACS) on board HST. [Larsen et al.[ ([2008[) constructed 
a CMD for the cluster NGC 1569-B and identified about 60 red 
supergiant stars (RSGs). This demonstrates that photometry can 
indeed be obtained for individual stars in clusters well beyond 
the Local Group, and motivated us to search for additional can- 
didates suitable for such analysis. Several other galaxies located 
within distances of a few Mpc have ACS imaging, obtained as 
part of our own and other programmes. We have inspected these 
images and looked for clusters that appeared sufficiently well re- 
solved into individual stars that photometry would be feasible. 
In this paper we present the resulting CMDs and compare with 
theoretical isochrones. Our aim is to quantify whether or not the 
observations can be well approximated by single isochrones. 

2. Data reduction and analysis 

2. 1 . Selection of the sample 

Even with the spatial resolution offered by HST, we cannot re- 
solve individual stars in clusters more distant than a few Mpc. 
For this work, we inspected images taken as part of our Cycle 12 
programme (GO-9774, P.I. S. S. Larsen). We found potentially 
suitable clusters in three out of five galaxies observed in this 
programme (NGC 1313, NGC 5236 and NGC 7793), while no 
suitable clusters were found in the two other galaxies (NGC 45, 
NGC 4395). The images obtained under this programme were all 
taken with the wide field channel (WFC) of the ACS, which has a 
pixel scale of 0'.'050 pixeL' . At the typical distances of ~ 4 Mpc, 
this corresponds to about 1 pc pixeL' so it is clear that we can- 
not obtain photometry in the cluster cores but only in the outer 
regions. The filters used for the ACS/WFC observations were 
F435W, F555W and F814W, roughly equivalent to the Johnson- 
Cousins BVI filters. For a detailed description of the data, see 
Moraetal. (2009). In addition, we include images of the clus- 
ters NGC 1569-A, NGC 1569-B and NGC 1705-1, taken with 
the ACS/HRC in the filters F330W, F555W and F814W as part 
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Table 1. Data for t he clusters analy se d in this paper. Refere nc es for distances: ' iMendez et al.l (l2002h . 
^Thim et"al ] (l2003h . ICai-ignanl (119851) . iGrocholski et al.l(l2008h . lTosi et al.l(l200lh 



Cluster 




R.A. & Decl. (2000.0) 


D [Mpc] 


V 


Ae 


log age 


Mass [Mq] 




FWHM 


1 


NGC1569-A 




04:30:48.08 


+64:50:57.2 


3.4" 


14.8 


2.30 


6.7 


6.6 X 10" 


~ 25 


a'057/0.9 pc 


1.12 


NGC1569-B 




04:30:48.88 


+64:50:51.2 


3.4" 


15.6 


2.30 


7.2 


1.2 X 10*^ 


~ 12.0 


a'084/1.4 pc 


1.04 


NGC1705-1 




04:54:13.48 


-53:21:39.4 


5.V 


14.8 


0.04 


7.1 


9.2 X 10' 


~ 12.5 


a'036/0.9 pc 


1.0 


NGC1313-F3- 


1 


03:17:47.76 


-66:30:18.7 


4.1° 


17.1 


0.47 


7.74 


2.8 X 10' 


5.9 


a'19/3.8pc 


0.83 


NGC5236-F1- 


1 


13:37:01.37 


-29:50:49.2 


4.5* 


17.1 


0.28 


7.44 


2.1 X 10' 


7.9 


ff.' 14/3.1 pc 


1.09 


NGC5236-F1-3 


13:37:04.67 


-29:50:35.2 


4.5* 


18.1 


0.28 


7.45 


8.1 X 10* 


7.9 


ff.' 15/3.3 pc 


1.13 


NGC7793-F1- 


1 


23:57:35.88 


-32:35:40.9 


3.3'- 


17.4 


0.07 


7.69 


9.2 X 10"* 


6.3 


a.'23/3.7 pc 


1.49 



NGC1313-F3-1 NGC5236-F1-1 NGC5236-F1-3 NGC7793-F1-1 



1000 




100 



6.5 7.0 7.5 8.0 8.5 9.0 

log t/yr 



Fig. 1. Number of stars with log g < 3.0 for an initial total stellar 
mass of IO^Mq. 



of programme GTO-9300 (P.I.: H. Ford). The F330W data are 
less suitable for stellar photometry due to lower S/N (especially 
for the heavily reddened NGC 1569), and we only make use of 
the F555W and F814W photometry in this paper We will from 
now on refer to the ACS F330W, F435W, F555W and F814W 
bands as U, B, V and /, although we do not actually transform 
the photometry to the standard Johnson-Cousins system. 

Table [T] lists basic data for the clusters we have exam- 
ined in detail. Coordinates were measured on the ACS im- 
ages; we caution that the absolute HST astrometry is only ac- 
curate to l"-2". Ages for the clusters in NGC 1313, NGC 5236 
and NGC 7793 were estimated from integrated (ground-based) 
UBV photometry (Larsen 1999,). using the e mpiri cal calibra- 
tion of UBV colour vs. age in iGirardi et al ] (11 99.51) which we 
note is strictly valid only for LMC-like metallicity. The clus- 
ters in NGC 1569 and NGC 1705 are significantly younger than 
those in the spirals, and age estimates from integrated light are 
more uncertain for these clusters due to the rapid, metallicity- 
dependent changes in integrated colours at very young ages. For 
NGC 1705-1 we initially adopt an age of 12 Myr, based on STIS 
spectroscopy (Vazquez et al. 2004). For NGC 1569-A we as- 
sume an age of 5 Myr (Ori gha et al.ll2001l: IWestmoauette et al.l 
20071) and for NGC 156 9-B we assume 15 Mvr dAnders et all 
2004 : iLarsen et al]|2008l) . The cluster masses were determined 
as a byproduct of the synthetic cluster experiments described be- 
low (Section [2.3b . We also list the main sequence turn-off mass 
(Mto) corresponding to the estimate d age of each cluster , as in- 
dicated in the Padua isochrone tables dBertelli et al.ll2009h . Since 
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Fig. 3. Similar to Fig.|2] but showing (B-V), My for those clusters 
where F435W band data are available. 



these refer to the bluest point reached during the core H-burning 
phase, the most massive stars still present in a cluster will be 
somewhat more massive. 



2.2. Photometry 

PSF-fitting photometry was done with the 
DAOPHOT/ALLFRAME package (IStetsonI 11994 . First, 
the FIND task in DAOPHOT was run on 1000 x 1000 pixels 
sections of the F555W images, centered on the clusters, and 
aperture photometry was obtained in all three bands with the 
PHOT task. Point-spread functions (PSFs) were constructed for 
each band based on typically 10-20 stars, and ALLFRAME 
was then run on all three images simultaneously to produce 
PSF-fitting magnitudes. In a second iteration of this process, im- 
proved PSFs were generated by subtracting all stars except the 
PSF stars from the images and redetermining the PSFs, and the 
FIND task was applied once again on the star-subtracted F555W 
images. Finally, ALLFRAME was run with the improved PSFs 
and combined object lists from the first and second iterations 
as input. The photometry was calibrated by obtaining aperture 
photometry of the PSF stars on the cleaned images and applying 
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Fig. 2. Top row: observed (V-I),Mv colour-magnitude diagrams for the seven clusters, showing all stars for which photometry 
was obtained. The red dashed curves indicate the approximate 50% completeness limits, while error bars show typical photometric 
errors. Second row from top: Colour-magnitude diagra ms for background a nnulus. Third row: Background-subtracted CMDs (see 
text for details). Bottom row: simulated CMDs, using iBertelli et all (l2009h isochrones for the ages and metallicities listed in the 
legends. 



(small) offsets to the PSF magnitudes to make them match 
the aperture photometry. The photometry was calibrated to the 
VEGAMAG system using zero-points published on the ACS 
web pages. We corre cted for foreground extinction according to 
ISchlegel et aU (T998), except for NGC 1569 where we assumed 
Ab = 2.30 mag (Israel 1988). 

Photometry was extracted in a circular region with a radius 
of 2" around each cluster, excluding the region within 5 pixels 
of the centre (0'.'25 for the WFC and 0'.'125 for the HRC) where 
crowding gets too severe to obtain useful photometry. Given that 
we are generally measuring stars located in the outskirts of the 
clusters, we expect some contamination of the cluster colour- 
magnitude diagrams by field stars. In order to correct for this, 
we used a large comparison annulus with an inner radius of typ- 
ically 7'.'5 and with a similar area as the region where cluster 
CMDs were obtained. For each star in the comparison annulus, 
the star in the CMD region which had the most similar B, V and / 
magnitudes was removed (allowing a maximum mismatch of 0.5 
mag). For clusters where no B data were available, the matching 
was done using V and / only. 

The top row in Fig. |2] shows the (V - I, V) CMDs for all 
stars within the central apertures, while CMDs for the compar- 
ison annuli are shown in the second row from the top and the 
decontaminated cluster CMDs are in the third row. Also shown 
in the third row are theoretical Padua isochrones for the ages 
given in Table [T] and metallicities indicated in each panel (see 



references in Sect. l2.4l l. Equivalent plots are shown in Fig.[3]for 
the (B - V, V) photometry. The cleaned CMDs all show a clear 
excess above the background and one can readily identify the 
red supergiant stars at V- / « 2 in several clusters (NGC1569-B, 
NGC1705- 1 , NGC 1 3 1 3-F3- 1 , NGC5236-F1 - 1 , NGC5236-F1 -3, 
NGC7793-F1-1). The CMD for NGC1569-A also shows a clear 
excess above the background, although this cluster appears too 
young to host a significant number of red supergiants. In the 
remainder of this paper we concentrate on the {V - I, V) pho- 
tometry, which constitutes a homogeneous dataset for the whole 
sample. 

A first rough indication of the 50% completeness limits (red 
dashed lines in Fig. |2]i and typical photometric errors was ob- 
tained by adding artificial stars to the cluster images, randomly 
distributed within the central aperture, and redoing the photom- 
etry. This was done for V magnitudes between 21 and 27 and 
V-I colours between -1 and -1-3, in steps of 1 magnitude and 
for 100 stars at each step. We required a minimum separation of 
5 pixels between two artificial stars. The completeness fraction 
was then determined as the fraction of added stars for which a 
star was recovered within a separation of 1 pixel. Since these ar- 
tificial stars are added on top of existing images, there is some 
probability of recovering a star even where none was added. This 
number was determined by looking for counterparts to the artifi- 
cial stars in the original input images and was typically on order 
~ 10%. An alternative approach would be to simply compare the 
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Fig. 4. Hess diagrams for Padua isochrones (iBertelli et alj|2009h 

of various ages and metallicites. The dotted line in each panel 
connects each point on the isochrones, while the shaded areas 
indicate the density of stars (assuming a Kroupa IMF). 

number of artificial stars added to an image with the additional 
number of detections. The problem with this approach is that a 
pre-existing star might be "masked" by an artificial star, leading 
to an underestimate of the number of recovered stars. It is clear 
that the classical artificial-star method is less than ideal for our 
situation and below we will adopt a more realistic approach to 
deal with the severe crowding in our fields. The classical artifi- 
cial star approach was also applied to the background annuli (see 
Fig. |2]i. These are less affected by incompleteness than the cen- 
tral apertures. Given the approximative nature of these complete- 
ness determinations, we have not attempted to take the differ- 
ent degrees of incompleteness into account when correcting for 
the background contribution, so we are likely over-subtracting 
it. However, from Fig. |2] and [3] it is clear that this mainly affects 
the CMDs near the completeness limits where errors are large 
anyway. The post-main sequence stars are, in most cases, hardly 
affected by the background correction, and well above the 50% 
completeness limits. 

2.3. Synthetic cluster experiments 

A simple comparison with isochrones is of limited use for in- 
terpreting data such as those in Fig. |2] since this fails to take 
into account the strongly varying density of stars along a par- 
ticular isochrone. More information can be gained by plot- 
ting the isochrones as Hess diagrams, as shown in Fig. |4] (for 
the Padua isochrones) and Fig. |5] (for the Geneva isochrones; 
[Leieune & Schaerer 2001). It should be noted that the Geneva 
isochrones are shown for the HST Wide Field / Planetary 
Camera 2 F555W and F814W filters, which are not com- 



Fig. 5. Same as Fig. ID but for Geneva isochrones 
dLeieune & Schaereij l200lh . Note that the Geneva isochrones 
are shown for WFPC2 filter definitions. 



pletely identical to the filters with the same names on the 
ACS. However, for the purpose of illustration the differences are 
marginal. 

The main characteristics of the Padua and Geneva isochrones 
are similar, although there are differences in the details. The pre- 
dicted Blue Hertzsprung Gap is evident at all ages relevant to 
the data analysed here. It is already clear from Fig. |2] that this 
feature is not nearly as evident in the data. In what follows, we 
will restrict the comparison to the Padua isochrones, for which 
photometry in the ACS filters is directly available. 

Even though the outskirts of the clusters appear resolved into 
individual stars in the ACS images, crowding still represents a 
major challenge and we cannot directly compare our observed 
CMDs with the Hess diagrams in Fig.|4]and|5] A realistic assess- 
ment of photometric errors and detection limits is crucial in or- 
der to carry out a meaningful comparison between the observed 
CMDs and model predictions. As we have discussed above, the 
classical approach of adding stars at random positions within 
an image is insufficient to fully capture the complicated effects 
of crowding in our data. In addition to the effects mentioned in 
Sec. 12.21 the background and crowding vary strongly with dis- 
tance from the cluster centre, and a different approach is needed. 
Rather than relying on "typical" completeness limits and pho- 
tometric errors, we therefore opted to generate artificial coun- 
terparts of the real cluster images, and carry out photometry on 
those. In this way, we could directly compare the CMDs of our 
artificial clusters with those of the actual clusters. 

To generate the artificial clusters, we first derived structural 
parameters for each real cluster using the ISHAPE task in the 
BAOLAB package (Larsen 1999). ISHAPE assumes an analytic 
model of each cluster, where for this work we use profiles of the 
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Fig. 6. Top panels: F555W images of the star clusters. Bottom panels: simulated clusters with the same apparent magnitudes, ages 
and sizes. Each panel measures 4" x 4". 



following type, found bv lElson et alJ (Il987h to give a good fit to 
clusters in the Large Magellanic Cloud: 



I(r) = k 



(1) 



Here, rc is a core radius and rj determines how steeply the en- 
velope falls off at large radii. ISHAPE convolves the analytic 
cluster model with the PSF and iteratively adjusts the structural 
parameters of the cluster model until the best match to the data is 
achieved. In our case, the signal-to-noise was sufficiently high to 
allow ISHAPE to solve for both and rj. The last two columns 
in Table □ list the best-fitting r] and FWHM (Full Width at Half 
Maximum) for each cluster, the latter related to rc and rj as 



FWHM = 2rc il^/i-^ (2) 

for the Elson et aP (Il987h profiles. For the clusters located in 
spiral galaxies (NGC 1313, NGC 5236, NGC 7793) the FWHM 
values are fairly typical com pared to other clusters of similar 
masses in spirals (lLarsenl2004l) . while the clusters in NGC 1569 
and NGC 1705 have extraordinarily compact cores. This is par- 
ticularly remarkable given that the latter are also the most mas- 
sive clusters in our sample, implying very high stellar densities 
in the centre of these clusters lO^M© pc"^). We do not tabu- 
late half-light radii, which are very uncertain because the 77 val- 
ues are close to unity for most of the clusters (the half-light radii 
are, in fact, undefined for rj < I, unless the profiles are truncated 
at some finite outer radius). 

Having determined the intrinsic intensity profile of each 
cluster, we then proceeded to sample stars at random from the 
model profiles. The profiles in Eq. ([T]l extend to infinity, but we 
imposed a cut-off at a radius of 50 pixels (2'.'5 for the WFC im- 
ages and 1'.'25 for the HRC images). Each star was assigned a 
mass drawn at random from a Kroupa (2002) IMF, and corre- 
sponding magnitudes in each band were looked up in the Padua 
isochrones. Isochrones were selected according to the cluster 
ages given in Table [T| and the total number of stars in each ar- 
tificial cluster was tuned such that the integrated V-band mag- 
nitudes matched those of the real clusters. The stars belonging 
to these artificial clusters were then added, one by one, to the 
actual images using the MKSYNTH task in BAOLAB. To the 
extent possible, we added the artificial clusters at locations that 
resembled those of the real clusters in terms of the background 
level and stellar density. The photometric analysis described in 
Sec. I2.2l was then repeated on the images with the artificial clus- 
ters added. 



Real and simulated F555W images of each cluster are shown 
in Fig. |6l Although our simulated clusters closely resemble the 
real clusters, a couple of caveats should be noted: First, while we 
have taken great care to model the photometric (in)completeness 
and errors as realistically as possible using simulated images, 
there are always simplifications involved in the generation of 
artificial datasets. For example, the PSFs used in the artificial 
cluster images are approximations to the actual PSFs of the par- 
ent images, so that stars in the artificial clusters will have slightly 
different PSFs than the real stars present in the images. This may 
introduce subtle differences in the goodness of the ALLFRAME 
fits, possibly affecting the photometric errors. This would be a 
concern especially if we used the same PSFs to generate artificial 
clusters and carry out PSF-fitting photometry for them. To alle- 
viate this concern and mimic the process by which PSFs were 
derived for real clusters as closely as possible, we re-derived 
the PSFs for the synthetic images using a set of synthetic ref- 
erence stars. These were placed in similar locations within the 
images as the real PSF stars in terms of crowding and back- 
ground level. Another potential concern is our implicit neglect 
of mass segregation. This effect causes massive stars to move 
towards the centre of a cluster and low-mass stars to move to 
the outskirts. However, since the light profiles of both the artifi- 
cial and real clusters are dominated by the massive (luminous) 
stars, we do not expect this simplification to affect our analy- 
sis very strongly. Potentially more problematic is our neglect of 
binary stars. While random blends are automatically accounted 
for, physical binaries are not. This is a complex problem: a full 
account of binaries needs to consider not just the blending of 
light from two stars, but also possible effects of stellar evolution 
in binary systems (mass transfer, etc.) for realistic distributions 
of mass ratios and orbital periods. We will comment further on 
some possible effects of binary evolution on the CMDs below 
(Sect. [331). 

2.4. Simulated versus observed CMDs 

CMDs for the simulated clusters are shown in the bottom rows 
of Fig. |2] and Fig. [3] Background contamination has been sta- 
tistically subtracted in the same way as for the real clusters. 
Also shown are the isochrones used to generate the lists of in- 
put magnitudes. We have somewhat arbitrarily, but guided by 
the CMDs in Fig. |2] made a division between "red" and "blue" 
stars at V-I = 0.8 In Fig.|2]we compare the observed and model 
luminosity functions for blue and red stars. The "turn-overs" in 
both the simulated and observed luminosity functions generally 
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agree fairly well with the approximate 50% completeness limits 
indicated in Fig. |2] 

For the artificial clusters, both the LFs and the synthetic 
CMDs are averages of 100 individual realizations of each cluster 
For each realization, a different random number seed was used 
for the sampling of the IMF and the cluster profile, and the centre 
of the cluster was shifted around randomly within a box of +50 
pixels. This approach allows us to put one-sigma error bars on 
the simulated LFs that include not only stochastic fluctuations in 
the clusters themselves, but also background variations. It should 
once again be emphasized that the luminosity functions have not 
been corrected for completeness effects, so we do not attempt to 
derive quantitative estimates of the shapes of these distributions. 
Instead, we keep the comparison with the simulated data strictly 
differential. Similarly, the simulated CMDs shown in the figures 
are "average" CMDs generated by first combining the 100 indi- 
vidual CMDs for each artificial cluster, and then randomly sam- 
pling 1% of the stars. 

Even though we employ PSF-fitting photometry, there is still 
the possibility that chance superpositions of two or several stars 
(as opposed to physical binaries) overlap so closely that they will 
be fitted as a single star. This is likely responsible for some of 
the data points in the synthetic CMDs that fall significantly off 
the input isochrones. The effect will be present both in the simu- 
lated and real data, and should not affect the relative comparison. 
Using our synthetic clusters, we can obtain some rough estimates 
of the amount of blending. The percentage of stars affected by 
blending will clearly depend on the specific criterion for a blend 
to occur If we (somewhat arbitrarily) define a "blend" as a star 
with a companion that is closer than 1 pixel and at most one 
magnitude fainter than the primary star, we find that roughly 
20% of the brightest stars in our CMDs are affected by blend- 
ing. However, this fraction depends very strongly on position 
and magnitude. 

In terms of the magnitude limits and scatter in the colour 
of main sequence stars, the simulated CMDs resemble the ob- 
served CMDs quite well, giving us some confidence that the 
modelling procedure captures the most important observational 
effects. However, a detailed comparison shows several signifi- 
cant differences between the modelled and observed CMDs. We 
will now comment further on these. We first discuss similarities 
and differences between the observed and model CMDs for each 
cluster, and then consider how various modifications to our ini- 
tial model assumptions might affect the model CMDs. 

NGC 1 569-A: Comparison with the other CMDs clearly shows 
this cluster to be the youngest in our sample, in agreement 
with analyses based on integrated light. It is actually a bi- 
nary cluster, as already hinted at in pre-ref urbishment mis- 
sion H ST Wide Field / Planetary Camera data (lO'Connell et al.l 
11994 . In our model CMDs we have assumed a metallicity 
of Z = 0.008 for this clus ter and NGC 1569-B (see dis- 
cussion in iLarsen et al.l l2008l) . Spectroscopic studies have re- 
vealed the presen ce of both Wolf-Ray et and red supergiant stars 
in NGC 1569-A ('OrigHa et al."2001). However, while the raw 
CMD in the top panel of Fig. |2]does show a number of RSGs, 
the background annulus shows a similar population so we find 
no statistically significant population of RSGs in NGC 1569-A 
itself. We do however note that a number of the brightest stars 
(at My ~ -9) appear somewhat redder than any stars appear- 
ing in the model CMD. These stars cannot be accounted for by 
any choice of isochrone: the isochrone used for the modelling 
(log(age)=6.65) is the oldest one for which no RSGs are present. 



For log(age)=6.70, the RSGs appear suddenly at V-I ~ 1.0. The 
observed and modelled LFs do not match particularly well ei- 
ther; on the blue side the simulated LF is significantly shallower 
than the observed one. The LF for red stars is most likely due to 
residual background contamination. 

NGC 1569-B: The observed CMD reveals a large number 
of RSGs (ab out 60) which have been discussed in detail in 
iLarsen et al.l f2008 ). Fig. Qshows that the shape of the luminos- 
ity function of the RSGs is fairly well matched by the simulation 
for the assumed age and metallicity, although the detailed mor- 
phologies of the simulated and observed CMDs do differ some- 
what. In particular, the observed luminosities of the RSGs show 
more scatter at fixed colour than the simulation. Overall, the sim- 
ulated CMD appears sparser than the observed one. Although, 
by construction, the total luminosities of the simulated clusters 
match those of the observed ones, the integrated luminosity of 
the measured stars is about 30% lower in the simulation for 
NGC 1569-B. Thus, the DAOPHOT photometry appears to be 
slightly less complete for the simulated cluster. This may be due 
to slight inaccuracies in the modelling of the cluster profiles or 
differences in the CMDs. 

Similarly to several other clusters, the simulated CMD shows 
a clear BHG which is not present in the real data. This difference 
is also clear when comparing the observed and simulated LFs. A 
small number of very luminous stars (My < -8) are observed; 
however, these are also present in the simulation indicating that 
they may well be blends or other artifacts due to the extreme 
crowding in this very dense cluster 

NGC 1705-1: The image of this cluster shows considerable 
substructure, with two "filaments" and another fainter cluster lo- 
cated about 1" from the centre of the main core. The metallicity 
of NGC 1705 is intermediate between those of the SMC and 
LMC (ILee & S killman 2004) and for the simulation we have as- 
sumed Z = 0.008. The observed CMD of this cluster is very 
similar to that of NGC 1569-B, indicating similar ages for the 
two clusters, and also for NGC 1705-1 the simulated CMD is a 
rather poor match to the observed one. Again, a BHG is clearly 
expected at My ~ -6, but not observed. Fig. [T] shows that the 
observed and simulated numbers of red supergiants match quite 
nicely, although the observations show the RSGs to be somewhat 
cooler and fainter than predicted by the models for this age and 
metallicity. 

NGC 1313-F3-1 : This is a remarkable object. With a core ra- 
dius of about 4 pc, this is a relatively extended cluster and the 
HST images show it to be resolved into numerous individual 
stars. The metallicity of NGC 1313 is 12-i-log O/H = 8.4 + 0.1 
([Walsh & Rov 1997 ), whi ch translates to [O/H] = -0.43 ± 0.1 
( Grevesse & Sauvall ll"998l) or Z = 0.007 assuming Solar abun- 
dance ratios. We thus use the Z - 0.008 Padua models for 
the simulated CMD, which to first order results in quite a good 
match to the data. This is one of the few cases where the ob- 
served CMD and LF show a hint of the BHG that is expected 
at My ~ -4. However, the observed CMD shows a spread of 
nearly 2 mag in the My magnitudes of the supergiant stars, a 
significantly larger spread than in the simulated CMDs. 

NGC5236-F1-1 and NGC 5236-F1-3: The metalhcity of most 
of the NGC 5236 disc is well above Solar dBresolin & KennicuttI 
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Fig. 7. Observed and simulated lumitiosity functions for red and blue stars in the clusters. 



and we have used the Z = 0.03 isochrones for the sim- 
ulations. However, the simulated CMDs would not look much 
different for the Solar metallicity isochrones (Z - 0.019). The 
integrated colours suggest identical ages for these two clusters, 
and their observed CMDs are indeed quite similar and in nei- 
ther case well reproduced by the simulations. NGC 5236-Fl-l 
is the brighter of the two by a magnitude, and the CMD shows 
about twice as many RSGs compared to NGC 5236-F1-3. The 
observed My magnitudes of the RSGs are about a magnitude 
fainter than the simulated ones, possibly an indication that the 
clusters are older than we have assumed. The mismatch between 
models and data is even worse on the blue side of the CMD. The 
models predict a clear BHG both in My magnitude and colour, 
whereas the data show no such separation. This problem cannot 
be solved by adopting an older age. 

NGC 7793-F1-1: The observed CMD is similar to that for 
NGC 1313-F3-1 but more sparsely populated, consistent with 
the ~ 3x lower estim ated total mass. We have assumed Z - 
0.008 jEdmunds & Pag el 1984). Many of the same remarks 
given for NGC 1313-F3-1 apply to this cluster, but at lower sta- 
tistical significance. Also here, there appears to be a significant 
1-2 mag spread in the My magnitudes of the supergiant stars. 

3. Discussion 

3.1. General comments on CMD comparisons 

Comparing our observed and simulated CMDs, we first note 
that they do agree in some respects. The colour of the main 
sequence is, in general, quite well reproduced, an indication 
that our adopted extinction values are reasonable. In addition, 
the CMDs show the red supergiants to be fairly well separated 
from the blue loop and main sequence stars, in agreement with 
model predictions for the relevant mass/age range (Fig. |4] and 
|5]l. However, there are also several areas of disagreement: 1) 
Even after accounting for observational errors, we confirm the 
presence of stars in the B HG where none are expected. This 
is a long-standing problem (IChiosilll998t ISalasnich et al.ll 19991) 



Table 2. Input for multiple-burst models, log f i - log indicate 
the logarithmic ages (in years) of the three bursts and /i - the 



relative amounts 


of stars 


formed 


in each 


burst. 






Cluster 


log ti 


log/2 


log t3 


/i 


fi 


/3 


NGC1569-B 


7.20 


6.95 


6.60 


0.70 


0.15 


0.15 


NGCl 705-1 


7.25 


7.05 


6.60 


0.45 


0.25 


0.30 


NGC1313-F3-1 


7.70 


7.50 


7.30 


0.65 


0.20 


0.15 


NGC5236-F1-1 


7.60 


7.20 


6.60 


0.85 


0.05 


0.10 



and has previously been not ed e.g. for the ~ 50 Myr old clus- 
ter NGC 1850 in the LMC dFischeret al.lll993h . Likewise, in- 
spection of the CMD for the 25 -Myr old cluster vdBO in M31 
(Perinaet al. 2009) shows no clear gap, although statistics are 
poorer here. 2) In several cases, the observations show a larger 
spread in the luminosities of the supergiant stars than expected 
from our artificial-cluster experiments. 3) The observed colours 
(i.e., effective temperatures) of the red supergiant stars generally 
tend to be somewhat redder than the predictions. This is seen 
both at super-solar metallicities (i.e., in NGC 5236), as well as 
in the more metal-poor (but younger) clust ers NGC 1569-B and 
NGC 1705-1 (see also iLarsen et alj|2008l) . The agreement be- 
tween models and data is best for the clusters NGC 1313-F3-1 
and NGC 7793-Fl-l, which are the oldest in our sample and 
probably have slightly sub-solar metallicities. 

In the following we discuss several modifications to our ini- 
tial, simple model assumptions that might affect the CMDs. 

3.2. Age spreads? 

In recent years, the traditional view of star clusters as simple 
stellar populations has been challenged by various pieces of ob- 
servational evidence. In addition to the long-standing puzzle of 
anomalo us abundan ce ratios of light elements in globular clus- 
ters (e.g. lKraftI 19791) . this includes multiple main sequence t urn- 
offs and subgiant branches in old GCs (Milon e et al.ll2010l) . as 
well as bi-modal or broad distributions of main sequence turn- 
off points in intermediate-age clusters in t he Magellanic Clouds 
dMackev et alJ l2008t iMiione et all l2009l) and the presence of 
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dual red clumps in such clusters (iGirardi et al.l2009l) . For ancient 
GCs there are few observational constraints on any actual age 
spread, and the observed CMD feat ures also appear to require 
strongly enhanced He abundanc es dNorrisI l2004t iPiotto et al] 
l2005HD'Antona & Venturall2Q07h . The LMC and SMC clusters 
where age spreads have been claimed all have ages of > 1 Gyr, 
and the implied age spreads (> 10^ years) would have dramatic 
effects on the CMDs of younger clusters. Theoretically, it is a 
challenge to understand how extended star formation could oc- 
cur within star clusters on long time scales. After a few Myr, the 
winds and radiation pressure from OB stars (or eventually, su- 
pernovae) are likely to evacuate any remaining gas, and the ma- 
terial subsequently returned by stellar evolution is insufficient to 
form the observed numbers of second-generation stars unless the 
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Fig. 10. Left panels: the division into brighter and fainter super- 
giant (SG) stars. Right panels: Spatial distribution of the fainter 
and brighter SGs. K-S tests are generally consistent with no dif- 
ferences in the radial distributions of the two populations, as in- 
dicated by the F-values in each panel. 



first generation has a very top-heavy IMF. This problem may be 
partly alleviated if a large fraction of the first generation of stars 
is preferentially lost via dynamical effects, and/or mass lo st from 
stars is supplemented by additional accreted material (e.g. BekkH 
1201 li) . It has been proposed that gas may indeed be accreted 
onto massive clusters as they pl ough through the interstellar 
medium of their parent galaxies dPflamm-Altenburg & Kroupal 
2009; Conroy & Spergel 2011]). 

So far, evidence of age spreads greater than a few Myr 
is scant in younger clusters. NGC 1850, one of the most lu- 
minous young clusters in the LMC, hosts a small number of 
stars with ages of only a few Myr that are located within a 
small subcluster (Fischer et al. 1993). Multiple stellar genera- 
tions (10-16 Myr and 30-100 Myr) have also been reported for 
the cluster Sandage 96 in the nearby Sc-type spiral NGC 2403 
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(IVinko et alJ l2009h although the older population in this case 
may be associated with the surrounding field. The giant Hii re- 
gion NGC 604 in M33 appe ars to host disti nct populations with 
ages of ~ 3 and ~ 12 Myr (lEldridge & Relafloll201 ll) . and star 
formation also seems to have occurred in the star forming region 
NGC 346 in the Small Magellanic Cloud over a per iod of 5-10 
Myr (iHennekemper et aT 2008; Cignoni et al. 2011) and possi- 
bly even up to ~ 20 Myr (De Marchi et al. 201 1). However, the 
latter two objects are probably too extended to evolve into single 
bound clusters. 

An age spread would result in diff'erent main sequence turn- 
off" points within a given cluster, thus providing a way to blur the 
CMDs, thereby potentially helping to fill the BHG and would 
also give a spread in the luminosities of the supergiant stars. 
We therefore now investigate to what extent this would help im- 
prove the agreement between our observations and the synthetic 
CMDs. 

As a starting point, we use the cluster NGC 1313-F3-1. As 
can be seen in Fig. |2] the observed colours of the supergiant 
stars are fairly well approximated by the models for this cluster 
It therefore appears plausible that a better overall match to the 
CMD could be obtained by includi ng a range of ages. We first 
use the CMD-fitting code FITSFH (" Silva-Villa & Larsenll20ldh 
to search for the combination of isochrones that provides the best 
match to the observed CMD. This code first constructs a Hess 
diagram for the observed data, and then attempts to match this 
by linearly combining Hess diagrams for individual isochrones 
(such as those in Fig. |4| after accounting for data incomplete- 
ness and observational errors. We used the errors estimated by 
ALLFRAME as input to FITSFH. Because of the large uncer- 
tainties in the completeness estimates, we have not actually in- 
cluded a detailed modelling of detection incompleteness across 
the CMDs in the SFH fits presented here, but simply restricted 
the fits to regions of the CMDs above the 50% completeness 
limits in Fig.|2l 

Fig. [8] shows the resulting star formation history (SFH) for 
NGC 1313-F3-1. The fit suggests a somewhat extended SFH, 
starting about 50 Myr ago and extending until 15-20 Myr ago. 
Integrating the SFH in Fig. [8] we find that about 2/3 of the stellar 
mass was formed at 7.6 < log t < 7.8, although most of this may 
have formed in a single burst at log f x 7.7. The remaining 1/3 
of the mass is in the tail extending to log t » 7.2. The first burst 
coincides quite well with the age derived from the integrated 
light. A similar analysis was done for three other clusters in our 
sample: NGC 1569-B, NGC 1705-1 and NGC 5236-Fl-l. For 
NGC 1569- A we do not have much freedom to experiment with 
different SFHs because the lack of RSGs essentially only allows 
the youngest isochrones, and we further excluded NGC 7793- 
Fl-1 and NGC 5236-F1-3 which have CMDs very similar to 
those of NGC 1313-F3-1 and NGC 5236-Fl-l, respectively, but 
with poorer statistics. 

While FITSFH will always output some best-fit star for- 
mation history, the quality of the "best fit" will be limited 
by the adequacy of the input assumptions, such as the stellar 
isochrones used, correct assessment of photometric errors, treat- 
ment of binaries, etc. For young stellar populations, it is gener- 
ally^ difficult t o obta in a^ery good fit to observed CMDs (e.g. 
[Siiva- Villa & Larsenll201 ih . and here we are faced with the ad- 
ditional difficulty of complicated incompleteness eff'ects and rel- 
atively small numbers of stars. We cannot, therefore, claim to 
put very tight constraints on the specific star formation histories 
of the clusters, but merely aim to explore whether an age spread 
might account at least in a qualitative way for some of the ob- 
served CMD features, and roughly how large it would have to be. 



Our final step is therefore to construct model CMDs based on the 
FITSFH output and compare these with the observed CMDs. As 
an approximation to the extended SFHs returned by FITSFH, we 
model each cluster as a superposition of three individual bursts 
with the ages and relative amounts of mass in each burst listed in 
Table |2l The relative weights of each burst were determined by 
integrating over the SFHs from FITSFH. We see that, in general, 
most of the mass is in the oldest burst. The implied age spreads 
are ~ 10-30 Myr Due to the presence of younger populations 
in each cluster, which decreases the overall mass-to-light ratio of 
the cluster, the total masses of the multiple-burst models are 25- 
40% less than those of the single-burst models listed in Table [T] 
The modelling as three distinct bursts was done mostly for com- 
putational convenience, as it allowed us to easily modify the set- 
up used for the single-burst models by simply adding three in- 
dividual "clusters" on top of each other. We could, of course, 
have increased the number of individual "bursts", but found that 
three individual sampling points in age provided an adequately 
smooth representation in the final CMDs. We emphasize that we 
do not intend to imply any preference for the clusters actually 
having experienced distinct bursts over smooth star formation 
histories. 

In Figure |9] we compare the multiple-burst models with the 
observations. The top row again shows the observed CMDs, 
while the middle row shows the CMDs for our multiple-burst 
models. Compared to the single-burst models, the fits are im- 
proved in the sense that we now get a larger range in luminosity 
for the supergiants, and the BHG is indeed blurred. However, 
the fits remain less than perfect. For NGC 1569-B, NGC 1705-1 
and NGC 5236-Fl-l, filling the BHG requires the inclusion of 
a younger isochrone (logfj - 6.60) that has no RSGs, so that 
the number of RSGs is decreased with respect to the single-burst 
models. At least for NGC 5236-Fl-l, this solution does not ap- 
pear physically plausible, and results in a number of spuriously 
bright blue stars in the model CMD. 

If the clusters indeed have an age spread, it would be of inter- 
est to study the spatial distributions of different sub-populations 
of stars within them. In Fig. [10] we attempt to make such a 
comparison, by dividing the supergiant stars into a brighter and 
fainter subsample in each cluster. The left-hand panels illustrate 
our adopted division for each of the four clusters, and in the 
right-hand panels we plot the spatial distributions of the stars. 
The circles indicate the areas in which photometry was carried 
out. In NGC 1569-B there is a hint that the brighter stars appear 
somewhat more centrally concentrated, but in the other clusters 
there are no statistically significant differences between the ra- 
dial distributions of bright and faint stars. This is quantified by 
the P-values from a K-S test, given in each panel. We do note that 
there is a curious asymmetry in the distribution of the brighter 
stars in NGC 1313-F3-1. The crossing time is only a few Myr, 
much less than the age of the cluster, so such substructure should 
be quickly erased by the internal motions of stars in the clus- 
ter If the brighter stars are as young as ~ 20 Myr, as suggested 
by our CMD modelling, and formed as a coherent group, this 
might help explain how some substructure could remain visible. 
Alternatively, the apparent asymmetry may simply be the result 
of a chance alignment. 

3.3. Extinction variations 

In the discussion so far, we have implicitly assumed that there is 
no significant spread in the extinction towards individual cluster 
stars. In the bottom row of Fig.|9]we show simulated CMDs that 
include a range of AAb -1.5 mag in the B-band extinction. As 
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in the multiple-burst models, this was done by simulating three 
distinct values {As - 0,0.75 and 1.5 mag). In these models, 
the ages of the model clusters are significantly younger than for 
the models in Fig.|2]and|3] By trial-and-error, we found that the 
best fits to the observed CMDs were achieved for input ages of 
logf = 7.0,7.0,7.50 and 7.25 for NGC 1569-B, NGC 1705-1, 
NGC 1313-F3-1 and NGC 5236-Fl-l, respectively, and for 70% 
of the stars belonging to the most extincted (Ab = 1.5 mag) bin. 

Although a spread in extinction will clearly contribute to 
blurring of the CMDs, it is evident from the figure that there 
are several difficulties with this scenario. In all clusters, most of 
the supergiants would have to be relatively strongly extincted, 
with only a smaller fraction of them scattering upwards in the 
CMDs towards lower extinction values. This, however, is in con- 
trast to the main sequence whose colours are well matched by 
isochrones with little or no additional extinction needed (see also 
Fig.|2]and[3]l. Furthermore, the small width of the main sequence 
in some clusters (especially NGC 1313-F3-1, and to some extent 
NGC 1569-B), also precludes a large scatter in the extinction 
values, unless the supergiants are systematically more subject to 
extinction variations than the main sequence. In conclusion, we 
find it unlikely that extinction variations have a major impact on 
the observed CMDs. 

3.4. Cepheids 

Another possible reason for the spread in the luminosities of 
the supergiant stars is that some of these stars are variable. In 
Fig. |9] we have indicated the location of the Cepheid instability 
strip (ISandage"etaDl2004 . While NGC 1569-B and NGC 1705- 
1 are probably too young to host Cepheids, it is possible that 
some stars in NGC 1313-F3-1 and NGC 5236-Fl-l (as well as 
NGC 7793-Fl-l and NGC 5236-F1-3) may be Cepheids. The 
pulsatio n amplitudes can am ount to more than 1 mag in the 
y-band (ISandaee et al.ll2004 . Although there is a hint that the 
spread in My magnitude is largest near the instability strip, a 
considerable scatter remains outside the boundaries of the strip, 
so it is questionable whether this can be the whole explanation. 
Observations of Cepheids might, however, shed light on possi- 
ble age spreads within the cluste rs via the period-age relation 
(lEfremovl200llBono et al.ll2005h . 

3.5. Binaries 

Binar ies are common among massive stars (e.g. iMason et al.l 
12009 ^ and can aff'ect the C MDs in various ways, see for example 
the excellent discussion bv lPols & MarinusI d 19941) . Although it 
is beyond the scope of this paper to model these effects in detail 
we will briefly discuss them. 

Unresolved multiples: The presence of an unresolved non- 
interacting companion can increase the brightness by up to a fac- 
tor two, i.e. change the magnitude by up to -0.75 mag. Although 
we must seriously consider the possibility that all stars have an 
unresolved companion, in most cases the companion will be con- 
siderably less luminous. Recent data suggest that the distribution 
of mass ratios (q) is roughly flat, implying that a typical compan- 
ion is only ha lf the mass of t he primary (Pinsonnea ult & Stanekl 
I2OO6 ; Kobul nickv & Frved2 007; Kiminki et al. 2009]). Given flie 
steep dependence of the luminosity on the mass, the effect on the 
brightness of such a companion would be less than -0. 1 mag for 
main sequence stars, and even less if the primary star has already 
evolved off the main sequence. If unresolved non-interacting bi- 



naries were to contribute significantly to any spread in the lumi- 
nosities of the supergiants, this would require that a significant 
fraction of these binaries have nearly identical masses so that 
both stars would be in the supergiant phase at the same time. 

For these clusters, where the stellar density is high, unre- 
solved higher order multiple systems may be rather common. 
These may either be physically bound triple or multiple systems 
or multiple stars that are aligned by chance, although the latter 
case is in principle taken into account by our artificial cluster ex- 
periments. However, even if the light of a star is polluted with 
the light of two additional equally bright stars the effect on the 
magnitude is still only -1.2 mag. This is comparable but still 
smaller than the size of the BHG as predicted by stellar mod- 
els, which is about ~ 2 mag. A shift of 2 mag is reached only 
for an unresolved multiple system out of which at least 6 stars 
are equally luminous. We conclude that unresolved multiple sys- 
tems may contribute to differences between our observed and 
modelled CMDs but can certainly not be the lone culprit. 

Interacting binaries: Data from nearby open clusters suggest 
about 4 out of 10 O stars have companions that are nearby 
enough to interact dSana & Evan s 2010). Even though the frac- 
tion of stars whose properties have been altered by binary inter- 
action may not be that large, the consequences of such interac- 
tions for the stellar luminosity are potentially severe due to the 
steep dependence of the luminosity on mass. In the most extreme 
cases, interaction can lead to a star that is twice as massive as the 
turn-off mass, which would imply an increase in the bolometric 
luminosity of a factor ~ 15 or ~ 3 mag. In practice, such cases 
should be rare and a more typical shift would be on the order of 
~ 5x, for a star that increases its mass by about 50% as a result 
of mass accretion or a merger of the original stars Q The effect 
on the CMD is still uncertain as it depends on various factors: 
the distribution of initial binary parameters, the amount of mass 
loss during interaction or mergers and the life time after interac- 
tion. In order to convert the changes in bolometric luminosity to 
magnitudes and colours, the effective temperatures of the merger 
products must also be known. 

Binary interaction also affects the post main sequence evolu- 
tion and therefore the luminosities of red and blue supergiants, as 
well as their relative numbers. The presence of a nearby compan- 
ion can strip the envelope from a star before it reaches giant di- 
mensions, effectively reducing the number of red supergiants by 
up to 40% (assuming the distribution functions of lSana & EvansI 
'2010). In addition, the accreting star or merger resulting from bi- 
nary interaction may under certain circumstances produce blue 
supergiants as suggested for the progenitor of SN1987A and the 
comp anion of SN1 993J jiHillebrandt et al1ll987l: IWalborn et all 
[19891 iMaund" et al.ll2004ft . This second effect depends, among 
other things, on the efficiency of internal mixing processes and 
the number and lifeti me of blue superg iants produced is there- 
fore hard to estimate dClaevs et al.ll201 1 1). However, in combina- 
tion the tw o effects may increas e the BSG/RSG ratio by up to a 
factor of 4 (Eldrid ge et al.ll2008h . 

In Figure [TT] we show how the CMD might be influenced 
if approximately two thirds of the stars are in interacting bina- 
ries. We use the binary population and spectral synthesis code 
BRASS (lEldridge & Stanwavll2009l) to create model stellar pop- 
ulations based on single stars (left) and binaries (right) at three 
different illustrative ages. The models have initial separations be- 
tween logia/Ro) = 1 and 4 in steps of 0.25 dex, while the mass 

' assuming that the luminosity of the star after interaction with the 
mass as L ~ ' 
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Fig. 11. Example CMDs for single stars (left panels) and binary 
(right panels) populations calculated with BPASS at three ages 
typical of clusters in our sample. Due to the nature of the bi- 
nary models the CMDs are created assuming a small age spread 
for the population of -I-/-0.05 dex. The grid cells are 0.1 by 0.2 
in y - / and My respectively. The grey scale is logarithmic in 
number of stars with the cell size. The dashed and solid con- 
tours represent when the density of stars in each cell is 0. 1 and 1 
respectively. 

ratio takes values of ^ = 0. 1, 0.3, 0.5, 0.7 and 0.9. We see that in 
general the effect of binaries, as discussed above, is to increase 
the luminosity of stars, including a population of "blue strag- 
glers" that populate the region above the MS turn-off, as well as 
a number of brighter blue and red supergiants. We note that the 
biggest increase is for the older populations shown here, where it 
produces an effect very similar to that of an age spread (Fig.|9]). 
These extra supergiants arise from binary mergers and secon- 
daries in binaries that become more massive by accreting mate- 
rial from their primaries. In both cases this leads to more mas- 
sive stars in the stellar population sometime after stars of such 
a mass would be expected to have exploded in supernovae for a 
simple single star population. This suggests that the creation of 



Table 3. Statistics of blue and red supergiant stars. The ratio 
M«sg/A^bsg is given for the observations, for the single-burst 
models (SSP) and for the multiple-burst models (Mult.). 



Cluster 


A'rsg 


A'bsg 




A'rsg 


/A'bsg 










Obs. 


SSP 


Mult. 


NGC1569-B 


69 


45 


1.53 


±0.09 


0.54 


0.28 


NGC1705-1 


29 


40 


0.73 


±0.03 


0.40 


0.16 


NGC1313-F3-1 


39 


32 


1.22 


±0.08 


1.42 


1.21 


NGC5236-F1-1 


20 


13 


1.54 


±0.30 


5.46 


2.79 



apparently younger stellar populations from binary evolution- 
ary paths might plausibly account for some of the differences 
between our observed CMDs and the models based on canoni- 
cal stellar isochrones (section [T2l l. While the code necessarily 
involves many simplifications in the modelling of binary evo- 
lution, and uncertainties in predicting the post-MS evolution of 
products of binary evolution are at least as large as for single 
stars, we note that the RSG to BSG ratios for these binary popu- 
lations in fact do provide values similar to those found in Section 
13.71 However, making more detailed predictions of the effects 
on the CMDs of these clusters is difficult due to the highly non- 
linear nature of binary stellar evolution, and beyond the scope of 
this paper 

3.6. Rotational mixing 

Massive stars are ob served with a wide range of rotation rates 
(iPennv & Giesll2009h . Rotation can induce mixing in the stellar 
interior during core hydrogen burning. This mixing brings fresh 
fuel towards the core of the star, which leads to an increased 
life time and to a l arger heli um core and thu s higher luminosity 
dMaeder & Mevne t 2000; B rottet all [20 Hal) . A distribution of 
initial rotational velocities will therefore lead to a dispersion of 
lifetimes and luminosities of stars of the same initial mass. 

However, the effect is expected to be rather limited. Even for 
very rapid rota tors, the life time is not increased by more than ~ 
30% dBrott et al. 201 lb), and less so for average rotators, which 
is smaller than the required age spread as discussed in Sect. 13.21 
The spread in post-main seq uence luminosity is predi cted to be 
similar, i.e., 0.3 magnitudes (iMaeder & Mevneil2000l) . which is 
again too small to account for the observed luminosity spread in 
our clusters. Furthermore, while rotational mixing can shift the 
red edge of the main sequence band slightly redwards, the same 
is tru e for the blue edge of the blue loops (iMaeder & Mevnet 
1200 Ih . such that the extent of the predicted blue Hertzsprung 
gap remains unaffected. 

3.7. Tlie red-to-blue supergiant ratio 

One long-standing problem in the modelling of massive stars is 
the number ratio of red to blue supergiants dLanger & MaedeJ 
1995). Determining this ratio observationally is far from 
straight-forward. Galactic open clusters generally contain only 
a small number of such stars, and in the field it is difficult to iso- 
late clean samples of red and blue supergiants because the latter 
are easily confused with massive main sequence stars. Due to 
the richly populated post main-sequence phases, the clusters in 
our sample offer a welcome opportunity to address this issue. 

In Table |3] we list the numbers of red and blue supergiants 
(Nrsg and A^bsg) found in the four richest clusters, using the 
colour division at (V - /) = 0.8 (as in Fig.|7]i and the magnitude 
cuts indicated by the lower dotted lines in Fig. [10] (in this com- 
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Fig. 12. Changes in integrated colours as a function of change 
in A^RSG /Nbsg ratio, relative to the ratio predicted by standard 
isochrones (A^rsg/A^bsg)std- 



parison we do not distinguish between the "bright" and "faint" 
subsamples). Specifically, the cuts are: 



NGC 1569-B 
NGC 1705-1 
NGC 1313-F3-1 
NGC 5236-Fl-l 



My < 0.5(y-/)-6.3 
My < 0.5(y-/)-6.1 
Mv < (V-I)-5.0 
Mv < (V-I)-5.Q 



Note that the samples of supergiant stars in Table [3] are larger 
than those found in any individual Galactic young clusters. We 
also list the number ratio A^rsg/A^bsg and its associated error, as- 
suming Poissonian statistics, along with the A^rsg/A^bsg ratio for 
our single- and multiple-burst simulations. The random errors on 
the simulated ratios are very small due to our averaging over 100 
realizations of each cluster The true uncertainty on the observed 
Nrsg/Nbsg is generally larger than indicated by the Poissonian 
errors, and typically dominated by the difficulty of isolating the 
BSGs due to the lack of a clear BHG. For the simulated CMDs in 
Fig.|2] the adopted cuts make a fairly clean selection of post-MS 
stars, but for the observed CMDs the number of BSGs depends 
on where exactly the magnitude cuts are made. 

For NGC 1313-F3-1 the agreement between the observed 
and simulated A^rsg/A^bsg ratio is excellent, but for the other 
clusters there is significant disagreement. For NGC 1569-B and 
NGC 1705-B, the observed ratios are significantly higher than 
the simulated ones, a discrepancy that is exacerbated in the 
multiple-age models. For NGC 5236-Fl-l it is very difficult to 
select a sample of BSGs in a meaningful way. 

3.8. Consequences for integrated colours 

Because the red and blue supergiant stars contribute a significant 
fraction of the light from a simple stellar population, it is of in- 
terest to consider how changes in this ratio will affect integrated 
properties such as broad-band colours. To assess this, we modi- 
fied the A^rsg /Nbsg ratios in the Padua models and re-computed 
various broad-band colours. We did this by adopting a division 
between red and blue supergiants at V - / = 0.8 and chang- 
ing the relative weights of stars on either side of this limit when 
computing integrated magnitudes. The (wavelength-dependent) 
integrated luminosity is then 



^Mmin 



w(M)LA{M)^{M)dM 



(3) 



whether a star is on the main sequence or is a red or blue su- 
pergiant. Here, ^ is the stellar IMF and Li(M) is the luminosity 
of a star with mass M. We kept the total number of supergiants 
constant, i.e.. 



/ 



wrsg I ^iM)dM + WBSG I ^(M)dM 

^{M)dM+ I ^(M)dM 



/ 

J Me 



which is the usual equation for a continuously sampled SSP apart 
from the factor w{M) which has different values depending on 



FigurefTSl shows the result of this exercise for two different ages 
(log? = 7.2 and logf - 7.75) and four different broad-band 
colours (U-B, B-V, V-I and V-K). This calculation was done 
for actual Johnson-Cousins colours. We plot the change in each 
colour with respect to the colour computed for the unmodified 
isochrones, as a function of the change in the A^rsg /Nbsg ratio 
(i.e., vvrsg/wbsg)- For B -V,V - 1 and V - K we see that the in- 
tegrated colours indeed do become redder when the A^rsg /Nbsg 
ratio is increased. However, for U-B the opposite happens. This 
somewhat counterintuitive result may be understood as follows: 
Taking the standard log f = 7.75 model as an example, the con- 
tribution of the RSGs to the integrated light is (0.8%, 4.0%, 13%, 
33%, 70%) in the (U, B, V, I, K) passbands, while the BSGs con- 
tribute (17%, 25%, 26%, 23%, 13%) in the same bands. Because 
both the blue and red supergiants are cooler than the main se- 
quence stars, the net effect of both is to produce redder integrated 
colours than would be the case for a cluster without any super- 
giant stars. At long wavelengths the red supergiants dominate the 
light, and increasing the Nrsq/Nbsg therefore has the expected 
effect of making the integrated colours redder However, at short 
wavelengths the contribution from the RSGs to the integrated 
light becomes very small, and only BSGs are able to produce a 
significant change in the integrated colours. Hence, increasing 
the Nrsg/Nbsg ratio in this case has the net effect of decreasing 
the overall contribution of supergiant stars to the integrated light 
at shorter wavelengths, thus making the integrated U-B colours 
bluer rather than redder 

Fortunately, the changes in integrated photometry appear rel- 
atively modest even for rather extreme (factor of 10) changes 
in the Nrsg/Nbsg ratio, at least for optical colours. U - B is 
the most age-sensitive of the colours shown in Fig. [12] (at least 
at young ages), and a change of A(U - B) of 0.1 mag corre- 
sponds to a difference of about 0. 1 dex in the age of a cluster It 
therefore appears that age estimates from integrated colours may 
not be too strongly affected by uncertainties in the RSG/BSG 
ratio. However, it would be worth investigating in more detail 
how other properties estimated from integrated colours (e.g. ex- 
tinction, metallicities, etc.) would be affected. In particular, the 
large effect on the near-IR colours might lead to substantial un- 
certainties in fits involving these colours that attempt to simul- 
taneously constrain age, metallicity and extinction. Combined 
with the fact that age determinations involving near-IR colours 
are also more sensitive to stochastic sampling of the stellar IMF 
(Fouesneau & Lancon 2010), this implies that the inclusion of 
more passbands (in particular extending to the near-IR) does 
not necessarily lead to a better determination of physical clus- 
ter properties from integrated colours. 

4. Summary 

We have presented resolved photometry for seven young mas- 
sive star clusters in five nearby galaxies, located at distances 
of 3-5 Mpc. While it is challenging to carry out photometry 
for individual stars in clusters at such distances, the richness of 
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these clusters allows us to obtain colour-magnitude diagrams for 
large samples of supergiant stars. In order to assess the effects 
of crowding, we have carried out extensive artificial-cluster tests 
where we generated synthetic clusters with properties as simi- 
lar as possible to those of our real clusters. By comparing the 
colour-magnitude diagrams of our real and simulated clusters, 
we reach the following conclusions: 

- We confirm previously known problems with canonical stel- 
lar isochrones: 

1. The standard models predict a clear separation in the 
CMD between the H-burning main sequence and the He- 
burning "blue loop" stars (i.e., the blue Hertzsprung gap) 
while no such gap is seen in the data. 

2. The observed red-to-blue supergiant ratio is in general 
not well reproduced by the models. However, the lack of 
a clear BHG makes it difficult to identify clean samples 
of blue supergiants for this comparison 

- We have investigated the effect of changing the red-to-blue 
supergiant ratio on integrated colours and find that even large 
changes in this ratio have relatively modest effects on opti- 
cal colours that would affect age estimates from integrated 
photometry by less than 0.1-0.15 dex. However, colours in- 
volving infrared passbands (such as V - K) are much more 
strongly affected, strengthening the conclusion that inclusion 
of such colours may actually increase, rather than decrease, 
the systematic errors w hen deriving cluster prop erties from 
integrated photometry ( Fouesneau & Lan9onll2010l) . 

- We observe a significant scatter (~ 2 mag) in the luminosities 
of the supergiant stars in some clusters. Comparison with 
our simulated clusters indicates that this scatter cannot be 
explained by observational errors alone. It may be partly due 
to variability, since some of the clusters are in the age range 
where Cepheid variables are expected. However, the scatter 
is seen also for stars outside the instability strip. 

- By including an age spread of 10-30 Myr in our model clus- 
ters we can achieve a better match to the observations. This 
is, however, much smaller than the spread that has been re- 
ported for older clusters in the Magellanic Clouds, and an 
age spread cannot fully explain the morphologies of all the 
CMDs. While uncertainties in stellar models may account 
for the differences between predictions and observations of 
the BHG, this is a less likely explanation for a spread in 
the supergiant luminosities. We discuss several other mecha- 
nisms that might contribute to differences between observed 
and model CMDs, and conclude that extinction variations 
and unresolved binaries are unlikely to be the main culprits, 
while interacting binaries might also have a significant effect 
on the CMDs. 

- In general, we find no significant differences in the spa- 
tial distributions of the brighter and fainter supergiant stars. 
In NGC 1569-B there is a hint that the brighter stars are 
somewhat more centrally concentrated, and in NGC 1313- 
F3-1 there is a curious asymmetry in the distribution of the 
brighter stars. 

Each of the clusters studied in this paper has its own 
unique characteristics and the current sample of seven objects 
is clearly still too small to draw general conclusions. However, 
our search of the HST archive is far from complete, and it is 
likely that several additional objects have already been observed 
that would lend themselves to similar analysis. Furthermore, 
the ongoing Panchromatic Hubble Andro meda Treasury sur- 
vey dDalcanton & PHAT Collaborationl201 li) will include a vast 



number of young clusters that may also be useful for addressing 
the issues discussed here (although most of them significantly 
less massive than those included here). Finally, it is our hope 
that the data presented here will prove useful as benchmarks for 
developers of stellar models. 
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